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Abstract 

A new asymptotic expansion method is developed to separate the Wheeler- 
DeWitt equation into the time-dependent Schrodinger equation for a matter 
field and the Einstein-Hamilton-Jacobi equation for the gravitational field in- 
cluding the quantum back-reaction of the matter field. In particular, the 
nonadiabatic basis of the generalized invariant for the matter field Hamilto- 
nian separates the Wheeler-DeWitt equation completely in the asymptotic 
limit of nip approaching infinity. The higher order quantum corrections of 
the gravity to the matter field are found. The new asymptotic expansion 
method is valid throughout all regions of superspace compared with other 
expansion methods with a certain limited region of validity. We apply the 

new asymptotic expansion method to the minimal FRW universe. 
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I. INTRODUCTION 



Recently the quantum field theory in a curved spacetime for matter fields has been 
studied in the context of quantum cosmology. In quantum cosmology one derives from the 
Wheeler-DeWitt (WDW) equation, a tentative quantum gravity theory, the quantum field 
theory for the matter fields in the curved space which is the time- dependent Schrodinger 
equation, a Tomonaga-Schwinger functional equation. The time-dependent Schrodinger 
equation for matter fields has an advantage over the canonical quantization on the solution 
space of wave equation for the matter fields in that one can take into account the higher 
order quantum corrections of the gravity to the matter fields and the quantum back-reaction 
of the matter fields to the gravity. Banks [|l| showed that the semiclassical limit (WKB ap- 
proximation) of the WDW equation for the pure gravity was the Einstein-Hamilton- Jacobi 
(EHJ) equation for the gravity, which turned out equivalent to the classical Einstein equation 
0. In Ref. the induced gauge potential due to some cosmological mode was considered 
in semiclassical gravity. The semiclassical gravity was extended to quantum cosmological 
models for the gravity coupled to matter fields by including the expectation value, a back- 
reaction, of the matter fields to the EHJ equation In a quantum cosmological model 
for the gravity coupled to matter fields whose mass scale is much smaller than the Planck 
mass, there are not only the back-reaction of the matter fields to the EHJ equation but also 
the geometric phases to the quantum states of the matter fields [p Hl^] - 

In this paper we develop a new asymptotic expansion method which separates the WDW 
equation for the gravity coupled to a matter field into the time-dependent Schrodinger equa- 
tion for the matter field and the EHJ equation for the gravity including the quantum back- 
reaction of the matter field. When one expands with respect to some basis of quantum states 
for the matter field and introduces a cosmological time the WDW equation is equivalent to 
a matrix equation which involves the cosmological time-dependent elements. In particu- 
lar it is shown in the nonadiabatic basis of the eigenstates of the generalized invariant for 
the matter field Hamiltonian that the WDW equation is equivalent to the matrix equation 
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which consists of dominant diagonal elements and nondiagonal perturbation elements pro- 
portional to the asymptotic parameter l/m^. In other words we are able to separate the 
WDW equation into the time-dependent matrix equation for the matter field and the EHJ 
equation with the back-reaction of the matter field. We obtain the higher order quantum 
corrections of the gravity to the matter field as a power series of l/rrip. We also compare 
the new asymptotic expansion method with other related works. Finally, we apply the new 
asymptotic expansion method to the minimal FRW universe. 

The organization of this paper is as follows. In Sec. II we introduce the new asymptotic 
expansion method for the WDW equation. In Sec. Ill the higher order quantum corrections 
of the gravity to the matter field are obtained. In Sec. IV we compare the new asymptotic 
expansion method with other methods used to derive the time-dependent Schrodinger equa- 
tion for the matter field. Finally in Sec. V we apply the new asymptotic expansion method to 
the Friedmann-Robertson- Walker universe minimally coupled to a free massive scalar field. 
The quantum back-reaction of the matter field and the EHJ equation are found explicitly 
using the generalized invariant for a well-known time-dependent harmonic oscillator. 



II. NEW ASYMPTOTIC EXPANSION METHOD 

With the fundamental constants such as c the speed of hght, h the Planck constant, and 
G the gravitational constant, inserted explicitly, the WDW equation takes the form 



*(/i„,0)=O, (2.1) 

„p UllaOllb ■ J 

where = c^/327rG is the Planck mass squared and Gab is the symmetrized supermetric 
on the superspace of the three-geometry. There are two asymptotic parameters 1/^ — > oo 
and rUp oo for the WDW equation, which can be compared with just one asymptotic 
parameter 1/h ^ oo for the Schrodinger equation in quantum mechanics. The asymptotic 
limit iTip — > oo can be considered as the limiting case either of G — > or ^ oo or both 
of them. We shall consider here only the weak gravitational coupling limit G — > without 
any loss of generality. 



The method in this paper is neither to evolve the wave functions of the WDW equation 
from some Cauchy initial data |T^ nor to separate the matrix effective gravitational field 
Hamiltonian equation by expanding the wave functions with respect to some basis for the 
matter field Hamiltonian [jl2[. The key point of our method is to separate asymptotically 
the Schrodinger equation for the matter field by introducing a cosmological time via the 
gravitational action which in turn includes the quantum back-reaction of the matter field. 
This can be achieved by setting the wave function in the form 

^>{ha, <i>) = exp [\S{ha)) $(0, ha). (2.2) 

Once we use the asymptotic parameter 1/^ ^ cxd in the WDW equation just as for the 
WKB expansion in quantum mechanics, then we obtain the following intermediate equation 

rGabi 2ih— — ir—— + ^ T^-r; 



$(0,/iJ=O, (2.3) 



where S is the gravitational action to be determined consistently later on and (a, b) denotes 
symmetrization with respect to the indices a and b. Assuming that the prefactor of the 
wave function in Eq. ( |2.2| ) oscillates and is peaked around a classical trajectory, we may 
first introduce the cosmological time along the classical trajectory in this oscillatory region 
of the superspace by 

or TTlp 0/l(a Ohb) 

and then rewrite Eq. (|2.3|) as 



inL _ ^^(vr^, 0, ha) + 2myVh ^'^'>R{ha) - ^'^ab^^^ 

. h 6^S h^ 6^ 

zGab-Ti ri 1" T: o^a 



2mp 6ha6hb 2mp 6ha6hb 



$(0,/ia) = O. (2.5) 

It should be remarked that Eq. ( |2.5| ) can be interpreted as the exact time- dependent 
Schrodinger equation for the matter field including the gravitational quantum back-reaction. 



since if we neglect the last two terms in Eq. ( |2.5|) , which is valid in the limit l/rrip 0, it 
is nothing but the conventional time-dependent Schrodinger equation plJ 



lh — ^<j), ha) = Hmi^^, 0, /ia)$(0, ha) (2.6) 

or 

provided that the gravitational action satisfies the conventional EHJ equation: 

1 AQ(0)/fC(0) 

7r-2Gat-ij--ir- - 2myVh ^''>R{ha) = (2.7) 
2m^ oha ohh ^ 

The second to the last term, —i{h/2mp)Gab{S'^S/6haShb), not only is asymptotically small 
as a power of l/rrip but also violates the unitarity of matter field equation, so we shall not 
consider this term in this paper. 

Our method to be developed below differs from others in that instead of taking the con- 
ventional EHJ equation in Eq. { \2.T\ ) to obtain the conventional time-dependent Schrodinger 
equation in Eq. (1^) we treat the exact time- dependent Schrodinger equation in Eq. ( [2.5[ ), 



expand the exact quantum states of the matter field by some basis, and finally solve the 
matrix equation equivalent to Eq. ( p. 51) . For the sake of simplicity, we shall use the ket- 
bra vector notation for the quantum states of the matter field, and act the bra vector on 
the ket vector to denote the inner product. The exact quantum state in Eq. ( p.5|) for the 
matter field can be expanded by some basis |$fc(0, ha)) which constitutes a complete set of 
orthonormal vectors as 

$(0,/l,) =^Cfc(M|$fc(0,/la)), (2.8) 
k 

where Ck{ha) are coefficient functions of the gravitational field only. Both Ck{ha) and 
\^k{4>, ha)) are still unknowns which will be determined systematically later on. Then Eq. 
( p.5| ) is equivalent to the following matrix equation 

ih — Ckiha) + n^^\ha)Ck{ha) 

Or 

+ E [-^'kliha) + n^uliha) + ^^in(^a)) C„(/.,) = 0, (2.9) 

where 



where 



2mp 6ha 6hb 



or 
dhadhf, 

+ Gab e-, e-, ^kn " 2i(A( 



6ha6hb 



6hb) 



(2.10) 



{Aa)kn = i ($fc(0, ha)\ JJ- ha)) • 



(2.11) 



The task of solving the WDW equation ( |2.1|) is now reduced to that of solving the matrix 
equation ( p.9| ). However, in the asymptotic limit of l/rrip 0, we may neglect the last term 
of summation in Eq. ( p.9| ) and get the approximate equation 



^n^Ck{ha) + n^'^Ckiha) + E (-^kliha) + 4n(/ia)) Cn{ha) = 0. 



(2.12) 



On the other hand, if one introduces a generalized invariant obeying the following invariant 
equation 



i{ha),Hm{ha) 







(2.13) 



for the matter field Hamiltonian, then there is a well-known decoupling theorem [15| for the 
generalized invariant such that 



(2.14) 



for different eigenstates, k ^ n, oi the generalized invariant: 



i{ha) |$fc(</), ha)) = Xk \'^k{(t>, ha)) 



(2.15) 



Since all of the off-diagonal equations in Eq. (|2.12 ) vanish, we have only the diagonal 
equations left 
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lh^Cu{ha) + (fi(°n/^a) - 4fc (/^a) + 4fc (M) Ck{K) = 0. (2.16) 

Recollecting that the gravitational action S is still undetermined, we may take a simple 
ansatz of the form 

fi^°HM-^^2(M = 0, (2.17) 

and 

th-^CkiK) + ol^^lihMK) = 0. (2.18) 



Equation (|2.17| ) is nothing but the familiar EHJ equation with the quantum back-reaction 



of the matter field included 

^ 'Gab^^ - 2myVh (=^)i?(/i,) + ($fe(0, K)\HU^^, (P, K) |$,.(0, K)) = 0. (2.19) 



2mp 5ha 6hb 

It should be noted that the EHJ equation (|2.19| ) is an implicitly coupled nonlinear equation 
in which the kth eigenstate is defined by the generalized invariant (p.l3| ), which in turn is 
defined using the cosmological time ( |2.4| ) through the gravitational action S. This means 
that the gravitational action does depend on the mode-number k, i.e. S'(fc), and so does 
the cosmological time r(fc). The physical implication is that there are infinite number of 
gravitational action S(^k) , where k runs over all quantum numbers of the generalized invariant 
which satisfies Eq. ( |2.19| ). The cosmological time r(-fc) is defined along each gravitational 



action S'(fc). It was shown in Ref. |]T1| that there is a spectrum of infinite number of wave 
functions that depend on the modes. The gravitational action S(^k) corresponds to the wave 
function and the cosmological time is defined along the classical trajectory around the 
peak of wave function. Now, the wave function (|2.2| ) becomes 

^ik)iha, 0) = exp Q5(fc)(/ij) <I>(fc)(0, K) (2.20) 

and so does the quantum state of the matter field 

$(fc)(0, ha) = Cn{ha) $(fc)n(0, ha)) . (2.21) 



We find the solution to Eq. ( p.l8| ) 

cf{K) = 4exp Q J (]S(/i,)dr(fc)) , (2.22) 

where dk is a constant. It was imphcitly assumed that all the Eqs. ( p.9|) through (p.l9|) 
should depend on the mode-number (A;); we shall, however, work with the {k) mode only 
and drop the mode-number hereafter throughout this paper. 

III. HIGHER ORDER QUANTUM CORRECTIONS 

We now turn to the time-dependent Schrodinger equation ( p.5|) for the matter field with 



the gravitational action 5* and thereby the cosmological time r determined by Eq. (|2.19|) 



and Eq. ( |2.4| ), respectively. The time- dependent Schrodinger equation, the gravitational 
action and the cosmological time are all defined along the (fc)th mode, whose mode-number 
will be omitted. Its matrix equivalent in Eq. ( p.9| ) now satisfies 

^hj^C,{ha) + n^hha)cM + ^E4'i(^a)c„(/l.) = 0. (3.1) 

By substituting 

Ckiha) = exp Q J n^^,}{ha)dT^ Ck{ha) (3.2) 
into Eq. ( p.l|) we obtain the following equation 

in-Cu{K) + ^E^S(^a)5„(/^a) = 0, (3.3) 



where 



^^^{K) = exp (-^ J nfi{ha)dr^ Qf^iK) exp J Q^^^{ha)dT^ . (3.4) 



The simple ansatz (|2.18|) was a solution to the asymptotic approximate diagonal equation 



( p.l6| ), therefore should be regarded as the zeroth order solution, 

~cf\K) = dk, (3.5) 



to Eq. (|3.1j ) with the EHJ equation (|2.19|) still satisfied. We may find the solution for Eq. 
( p. 3D perturbatively in a power series of li/2m'^: 




where they satisfy the recursion equation 

^^~^f{ha) + E ^^kKKYct^Kha) = 0. (3.7) 

" n 

The solution to the recursion equation is the ordered /th multiple integral 

rii ,i=l,...,l 

For an instance, when an initial data c^^ = dk^kn, dk = constant, is imposed, the quantum 
state of the matter field including the gravitational correction up to the first order is 

\ p n / 

(3.9) 

Thus the first order transition to the other states from the initially prepared state is sup- 
pressed by a factor of l/rrip because there comes a factor of h from the integration and 
similarly the nth order transition is suppressed by a factor of (l/m^)". It is to be noted 
that the procedure employed in this paper is quite similar to the perturbation method for 
a quantum system with a small perturbation term; here the parameter for the smallness 
of perturbation term is h/2m^ and Eq. ( |3.8| ) is the perturbative series for the given choice 
of the basis of eigenstates of the generalized invariant, which is also the solution for the 
unperturbed equation (|2.6D. 



IV. COMPARISON WITH OTHER RELATED WORKS 

In this section, we shall compare the result of this paper with other related works. 
First, we shall answer in part the question whether the conventional EHJ equation ( p.7| ) 
or the EHJ equation (p^.l9D is right through the investigation of a quantum cosmological 



model. Suppose a quantum cosmological model with the gravitational super-Hamiltonian 
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Hg{K) = ^G'afevr^TT^ - 2ml^Vh ^^^R{K). (4.1) 

When coupled to the matter field, the gravitational super-Hamiltonian ( [4 .11 ) leads to the 
WDW equation Quantizing the gravitational Hamiltonian Eq. ( [4.1| ) by substituting 

Tx"" = {h/i){6/6ha), taking the wave function of the form \l/(/ia) ~ exp{iS{ha)/h), and keeping 
dominant terms only, we just obtain the conventional EH J equation 

^ Ga,^^-2myVh('^R{h^) = 0, (4.2) 



2m2 ShaShb 



which is the same as Eq. (|2.7|) . Finally, putting the conventional EHJ equation ( |4.2|) 
and neglecting the last two terms in Eq. (|2.5| ), we obtain the time-dependent Schrodinger 
equation 

Z/I — $(0, ha) = HmiTTcj,, 0, /ia)$(0, (4.3) 

or 

We found the exact quantum states in terms of the eigenstate of the generalized invariant 



H<P,ha) = E 4exp {^-jnfMdr') \^k{<P,K) > . (4.4) 
It is to be noted that the conventional EHJ equation ([4.2[) is equal to the approximation 



^^^\ha) = and the gravitational field-dependent coefficients in Eq. ( [4.4| ) can also be 
determined by directly integrating Eq. ( p.l6|) with both n^'^\ha) = and n^^\ha) = 
substituted. However, the lowest order solution (|2.18|) of the coefficient function does not 



depend on the gravitational field, i.e., is a constant and the higher order solution ( p.8|) 
depends on the gravitational field. The significant difference on the coefficient functions 
comes from the fact that in this paper we have used the EHJ equation (|2.19|) with the 



quantum back-reaction of the matter field rather than the conventional EHJ equation ( [4.21) . 



Below we put forth a criterion on whether Eq. ( |2.19| ) or Eq. ( |4.2| ) should be used for the 
correct gravitational action. 



In order to show the relation between the EHJ equation ( |2.19| ) and the conventional 



EHJ equation ( [4.2| ), we expand the action ( p^.l9D perturbatively in the inverse power of the 
Planck mass: 
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m. 



2{l-n) (^{n) 



k 1 



(4.5) 



whose lowest order action obeys 



- 2c^^^h ^^^R{ha) = 0, 



(4.6) 



and the first two higher order actions obey 

-,(0 

0(m°) : 



0(m-') : G 



ab 



^2 ' 



5h(a ^h) 



2 6hr, 6hh 



0. 



(4.7) 



It is worthy to note that the lowest order contribution (|4.6| ) gives nothing but the conven- 
tional EH J equation ([4.2|) . So the above question whether Eq. ( 2.19 ) or Eq. (12) should be 
used is closely related to the question whether the asymptotic expansion of the gravitational 
action (|3) gives the correct gravitational action. 

However, contrary to the belief widely accepted that the conventional EHJ equation ( [4.2| ) 
gives the correct gravitational action, we give a counter-example showing that the asymp- 
totic expansion ( ^4.5| ) leads to a wrong gravitational action. In the case of the Friedmann- 
Robert son- Walker universe minimally coupled to a scalar field with a power-law potential, 
the WDW equation takes the form 



2mp a da? 



from which it follows that Gn 



- 2mlc^a - h^ — — + 2a^U{(f)) 



-1/a and 



equation ( |4.2| ) becomes 



1 1 (dS^ 



2m2 a \da 



4- 2m^(?a 



0. 



0. 



(4.8) 



+ 2a^U{4). The EHJ 



(4.9) 



Direct integration by quadrature yields the gravitational action S{a) = -^irri^ca? and the 
wave function 



(4.10) 



\E'(a, 4) = exp { ±—mpCa | $(0, a 
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Likewise, the lowest order gravitational action (|4.6| ) also has the same value S^^^ = ±.ica?. 
Therefore, the cosmological time Eq. (|2.4| ) leads to an imaginary one r = ±iln(a/2c). Both 
the conventional EHJ equation (|4.2| ) and the dominant term of the asymptotic expansion of 
the gravitational action (|4.5|) always lead to the wave functions with an exponential behavior 
due to the curvature term, the second term in Eq. ( [4.8|) . However, both in the adiabatic 
basis method [|T^ which expands the wave functions by the gravitational field-dependent 
eigenfunctions of the matter field Hamiltonian and in the superadiabatic expansion method 



in which transitions among different eigenstates are taken into account during the evolu- 
tion of the Universe, the resulted wave functions show not only the exponential behavior for 
a large three-geometry but also the oscillatory behavior for an intermediate three-geometry 
depending on the quantum number of the matter field Hamiltonian. The oscillatory behav- 
ior of wave functions is inevitable to the classical Lorentzian universe such as the present 
Friedmann-Robertson- Walker universe. Therefore, in the case of the minimally coupled 
Friedmann-Robertson- Walker universe it can be inferred that the conventional EHJ equa- 



tion ( [4.2|) has a certain limited region of the large three-geometry for the validity, whereas 



the EHJ equation (|2.19|) holds not only for the large three-geometry prevailing with the 



curvature term but also for the intermediate three-geometry prevailing with the quantum 
back-reaction of the matter field. The EHJ equation ( |2.19| ) should be used in order to give 
the correct gravitational action valid for all the regions of superspace. 

Second, we shall compare the new asymptotic expansion method with the adiabatic 
expansion method. Quite similarly as in the new asymptotic expansion method in which 
one expands the quantum state of matter field by the nonadiabatic basis of the eigenstates of 
the generalized invariant, in the adiabatic expansion method one expands the quantum state 
of the matter field by the adiabatic basis of the instantaneous eigenstates of the matter field 
Hamiltonian itself and includes the quantum back-reaction of the matter field with respect 
to the adiabatic basis. By defining the instantaneous eigenstates 

Hm{ha) ha)),^ = Hk{ha) K)) , (4.11) 
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and expanding the quantum state 

(^{(j), ha) = ^ CadAha) |$fc(0, ha)) , (4.12) 
k 

one obtains the matrix equation 

ih — Cad,k{ha) + ^^°\ha)CadAhcL) 

Or 

+ E [-^alM + ^%Aha) + ^,n%^{ha)^ CadAha) = 0, (4.13) 



where 



n^'\ha) = -T^Ga,^^ + 2myVh (''>R{ha), 

Zrrip ohaOhb ^ 

^alkAK) = ad ($fc(0, ha)\HAT^<p, 0, ha) |$„(0, ha)) ad , 
^alkAha) = ih ad ($fc(0, Ml ^ |$n(0, 

or 

^%n{ha) = Ga, ad ($ fc(0, /^a) | |$n(0, ha)) ad 

kn '^^(^■^(ad,a) kn r-, ; 

(4.14) 

ohadhb dhb) 



ad ' 



where 

{Aad,a)kn = iad ('J'fc(<^, ha)\ — |$„(<^, /la))arf ■ (4.15) 

Again, in the asymptotic hmit of l/frip 0, we may neglect the last term in Eq. ( [4.13| ) and 
get the approximate equation 

th^CadAha) + n^'\ha)CadAha) + E {-^alkniha) + ^fl,kniha)) C^ha) = 0. (4.16) 

It should, however, be remarked that in the adiabatic expansion method the ofF-diagonal 
elements of the coupling matrix do not vanish, Q^f}^ ^ Q^^^ for A; 7^ n, in strong contrast with 
those in the nonadiabatic expansion method. Therefore, one should solve the whole adiabatic 
approximate matrix equation ( |4.1(j| ) with the elements of the coupling matrix accounting for 



transition between different eigenstates instead of the frequently used adiabatic diagonal 
equation in Eq. ( [4.16| ) 
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5t 



(4.17) 



Because f^^]^'*;.^ for k ^ n, have an order of magnitude comparable to ^ad,kk^ not justified 
to use the adiabatic diagonal equation ([4.17|) , which is frequently used in literatures under 
the assumption that the off-diagonal elements be neglected. From Eq. ( fl.l7| ) one also 
obtains the frequently used EHJ equation with the adiabatic quantum back-reaction of the 
matter field 



1 G^,^^-2myVh^'^R{K 



2mp 6ha 6hh 



ad{'^ki^,ha)\HmiT^^A,ha) |$fc(0, /la))„rf = 0. (4.18) 



Here we have the same subtlety as in Sec. 11 that the EHJ equation with the adiabatic 
quantum back-reaction depends on the mode-number k. So we have the gravitational action 
S(^k) and the cosmological time T(^k)- The difference between the EHJ equation ( [4.18| ) and the 
EHJ equation ( |2.19| ) is that in the former the quantum back-reaction of the matter field is 
explicitly given, whereas in the latter it is implicitly determined via the cosmological time 
which is defined by the gravitational action. 

Third, there has been a recent study of geometric phases as a mechanism for the asymme- 
try of the cosmological time . In particular, in the basis of eigenstates of the generalized 
invariant one may define a gauge potential (Berry connection) 

5 



dha 



>, ha), 



(4.19) 



where 



^ |<l>o(0, ha) > ^ 
\^l((j),ha) > 



|$n(0, ha) > 



(4.20) 



V ■ / 

is a column vector, and the asterisk * and the superscript ^ denote dual and transpose 
operations, respectively. One can show that 
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and 



^[7(0, ha) = -Z^Gab-^Al^U{(P, ha), (4.21) 
OT ?7ip Oll(^a 



Gab^4r^{<P, ha) = Gab f ^ + {AaA,)A [7(0, ha). (4.22) 

ohaOhb \6ha J 



Then the matrix equivalent ( |2.9| ) to the WDW equation is entirely determined by the grav- 
itational action and the gauge potential. It should be remarked again that the matter field 
Hamiltonian gives not only a back- reaction to the EH J equation (p. 19 ) but also the geometric 
phase term 

($,(0,Mh^^|$,(0,M) = ^Gab^(A,^ (4.23) 
OT dh(^a ^ ^^'^ 

to the quantum state of the matter field. It should also be noted that the gauge potential 
defined in terms of the eigenstates of the generalized invariant does always give nontrivial 
diagonal elements in strong contrast with the gauge potential defined in terms of instanta- 
neous eigenstates of the Hamiltonian whose diagonal elements vanish for real eigenstates. 
On the other hand, when the gravitational field Hamiltonian ( [4.1|) acts on the wave function 
expanded by the eigenstates of the generalized invariant as 

^{haA)=U^{(t>,ha)-^{ha), (4.24) 

it becomes the matrix nonadiabatic gravitational Hamiltonian equation and acquires the 
induced gauge potential ( [4. 19] ) 



Hg^iha) 



^Gab{r - Aa){7r' - Ak) - 2myVh '-'^Riha) 



2m2 



^{ha). (4.25) 



It was shown that there was a remarkable decoupling theorem as mentioned earlier canceling 
the off diagonal gauge potential and the expectation value of the matter field Hamiltonian 
at the classical level. After the matter field Hamiltonian is included, the total Hamiltonian 
Hg + Hm acting on the A;th wave function 

^U{ha, 0) = exp (^J^Sk{ha)^ |$fc(0, ha)) (4.26) 
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leads to 



2ml ShaShb 



(a 



2ml 



- 2myVh ^'^R{K) + ha)\Hm{^^, 0, |$fc(0, K)) = 0. (4.27) 



This equation can be rewritten as 

[Wb'"^') -2^^'^+('^fc('^'Ml^m(7r^,0,MI*fc(0,M) = O. (4.28) 

Eq. (f4.27|) is the diagonal equation of Eq. ( p.9| ) with the wave function ( |4.26| ) and the zeroth 
order solution Ck = constant. In our new asymptotic expansion of the WDW equation what 
corresponds to the gauge potential is resigned to the matrix equation ( |2.9| ) for the matter 
field. 



V. MINIMAL FRW UNIVERSE 



We consider in detail the Friedmann- Robert son- Walker universe minimally coupled to a 
free massive scalar field, whose WDW equation is 



2m c a — h 

2ml ^ ^ '^^ 



+ m^a^d)^ 



0. 



(5.1) 



Here m is the mass of the scalar field and the matter field Hamiltonian is given by 

1 



Hm = -tVT^ + m a 



2^3 j2 



(5.2) 



In order to find the generalized invariant, the first thing to do is to find the classical equation 
of motion for Eq. (|5.2|) 



{T)+3^<P{T)+4m^^{T) = 0, 
a(r) 



(5.3) 



where the cosmological time will be determined later on through Eq. ( |2.4| ). The cosmological 
scale factor a and the scalar field depend implicitly on the cosmological time. Under the 
assumption that the classical solutions, 0i(t) and 02 (t), to Eq. ( p. 3] ) are given explicitly, it 
is known that the generalized invariant is given by [|18| 
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Hr) = 9-{r)^ + 9o{r) ^ " + <7+(r)^, (5.4) 



where 



^_(r) = Ci0i(r) + C20i(r)02(r) + 03^2 (r), 

9o{r) = (ci0i(r)0i(r) + | (0i(r)02(r) + 0i(r)02(r)) + c,Ur)Mr)) , 

= (ci0?(r) + C20i(r)02(r) + C30^(r)) . (5.5) 

We may introduce the cosmological time-dependent creation and annihilation operators 
of the generahzed invariant 



where 



^0 = ^Jg+ir)g4r)-gi{T) (5.7) 
is a constant of motion. The generahzed invariant can be written as 

i{r)=ujo(b^{r)b{T) + ^y (5.8) 



The eigenstates are the number states 



\n,T) = ^b^-{r)\0,T), (5.9) 



where the ground state is annihilated by the cosmological time-dependent annihilation op- 
erator 

6(r)|0,r) = 0. (5.10) 



In the number state representation the gauge potential ( [4.19| ) reads that 

A{t) = za{r) (6t(r)6(r) + ^ + ^ (/3(r)6^(r) - P*{rp') , (5.11) 

17 



where 



a(r) 



1 g_{r) d 



(r) 



2i loq dr yg-ir) J ' 
1 1 dg^ir 



2g4T) dr 

The gauge potential ( |5.11| ) has the matrix notation 

1 



+ a r 



(5.12) 



n-2 



■,n+2i 



- -/3*(r)v/(n + l)(n + 2)4. 
and the gauge potential squared 

AUr) = ((-a^ + {n' + n) + ^PP* - ^a'^ 

- ^aP{2n - l)^n{n-l)6k,n~2 + \aP*{2n + 3) ^'(n + l)(n + 2)4,n+2 

- - l)(n - 2){n - 3)Sk,n-4 

+ ^rV(n + l){n + 2){n + 3)(n + 4)4, 

One can show that 

2K + ^2(^)) + 2mW(r) 



(5.13) 



:,n+4- 



(5.14) 



2ujQa^g(T) 
^ 9-{t) 9 / go{r) " 
2 uo dr \g-{T)^ 



2 



(5.15) 
(5.16) 



.OA 



kn 



a da? 



dr 



d f a 



A 



kn 



6kn - 2iA 



kn 



d_ 

da 



(5.17) 



fis are functions of the cosmological time (|2.4|) through the dependence of the classical 
solutions ( |5.3| ) and thereby the generalized invariant ( ^.5|) on a(r). Therefore Eq. (p.l9|) 
which now reads 

2 



' -(^] -2mic^a + ^ir) 



2mp a \ da 



0, 



(5.18) 



together with the definition of the cosmological time 
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determines the action S ClS db function of a(r). 



VI. SUMMARY AND DISCUSSION 

In summary, we have developed a new asymptotic expansion method according to 
which the Wheeler-DeWitt equation was separated into the Einstein-Hamilton- Jacobi equa- 
tion with the quantum back-reaction of the matter field included and the time-dependent 
Schrodinger equation for the matter field. In the new asymptotic expansion the Wheeler- 
DeWitt equation was equivalent to two coupled nonlinear functional equations consisting 
of the cosmological time (|2.4]) and of the time-dependent Schrodinger equation (|2.5| ) for 
the matter field or the matrix representation ( p.9| ) of it. The time-dependent Schrodinger 
equation or the matrix equation for the matter field included quantum gravitational correc- 
tions. In particular we have found the exact quantum state ( |2.8| ) in the nonadiabatic basis 
of eigenstates of the generalized invariant for the matter field Hamiltonian by solving the 
matrix equation whose solution consists of Eq. ( |2.17D and Eqs. ( p. 61 ) and (|3.8|). The zeroth 



order quantum state in Eq. (|2.8|) was the nonadiabatic basis itself, which motivated the use 



of the generalized invariant to solve the time-dependent Schrodinger equation in Ref . . It 
was found that the quantum corrections of gravity gave rise to transition of quantum states 
of the matter field whose first order transition rate ( p.9| ) was suppressed by the factor of 



l/m^, and the nth order transition rate by the factor of {l/mT)'^. Moreover, Eq. ( p.l7|) was 



nothing but the Einstein-Hamilton- Jacobi equation ( p.l9| ) with the quantum back-reaction 
of the matter field included. The Einstein-Hamilton- Jacobi equation is an implicitly coupled 
nonlinear equation for the gravitational action by which the cosmological time ( |2.4| ) used to 
define the generalized invariant is defined. Since the quantum back-reaction of the matter 
field depends on the mode-number, the cosmological time as well as the gravitational action 
do depend the mode-number. The physical meaning of the mode-dependent gravitational 
actions is that the time-dependent Schrodinger equation (p.5|) or its matrix equation ( p.9|) 

19 



should be defined along the corresponding gravitational action. In fact there are infinite 



number of wave functions which are peaked along the gravitational action [0]. The higher 
order quantum corrections of the gravity to the matter field in Sees. Ill and V are calculated 
along a specific mode-number gravitational action. 

It has been shown through an investigation of the minimally coupled Friedmann- 
Robert son- Walker universe that the Einstein-Hamilton- Jacobi equation ( p.l9|) gives indeed 
the correct gravitational action rather than the conventional Einstein-Hamilton- Jacobi equa- 
tion ([4.2|) providing an oscillatory regime necessary for the emergence of the Lorentzian 
universe. The new asymptotic expansion method based on the generalized invariant has the 
advantage that one has already canceled the off-diagonal terms between two different eigen- 
states due to the remarkable decoupling theorem of the generalized invariant, whereas one 
has to take care of the off-diagonal terms between two different instantaneous eigenstates 
of the matter field Hamiltonian in the conventional adiabatic expansion method. Further- 
more, by introducing the gauge potential (Berry connection) ( [4. 191 ) we were able to express 



explicitly the gravitational correction ( ^.221 ) as well as the back-reaction of the matter field. 



Finally, we have applied the new asymptotic expansion method to the Friedmann- 
Robertson- Walker universe with a minimal scalar field. The generalized invariant was found 
in terms of the classical solutions of the matter Hamiltonian, from which the back-reaction 
( ^.15| ), the geometric phase (|5.16| ), and the coupling matrix (|5.17|) of the gravitational cor- 



rections are derived explicitly. 

Considering an analogy, if any, between a quantum cosmological model of gravitational 
field and matter fields and a quantum system of heavy particles and light particles, the 
gravitational field and heavy particles behave as slow variables, and the matter fields and 
light particles behave as fast variables [|T^. In the quantum jargon of the fast and slow 
variables the new asymptotic expansion method provides us with a very systematic method 
to separate one equation for the fast variables and the other equation for the slow vari- 
ables. The fast variables obey a parameter-dependent quantum mechanical equation whose 
parameter is determined by the slow variables and the slow variables satisfy a classical 
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Einstein-Hamilton- Jacobi equation with the quantum back-reaction of the fast variables. In 
quantum cosmological models the gravitational field with the Planck mass scale behaves as a 
heavy particle obeying classical Einstein-Hamilton- Jacobi equation with the quantum back- 
reaction of the matter fields which consists of the gauge potential as well as the expectation 
value of the matter field Hamiltonian. The time- dependent Schrodinger equation for the 
matter field can also be solved using the generalized invariant. The free massive scalar fields 



in quantum cosmological models have been solved in Ref. [|I^ and have provided a mech- 
anism for the cosmological entropy production during an expansion and recollapse of the 
Universe. It is a new feature of the new asymptotic expansion of the Wheeler-DeWitt equa- 
tion that the geometric phases are an inevitable consequence to the quantum states of the 
matter field as well as the expectation value of the matter field Hamiltonian to the Einstein- 
Hamilton- Jacobi equation which is the counterpart of the classical Einstein-Hamilton- Jacobi 
equation. The gauge potential in the nonadiabatic basis can not be gauged away even for 
the real eigenstates. However, the gauge potential can always be gauged away for the adia- 
batic basis of real eigenstates of the matter field Hamiltonian in the conventional expansion 
method. Thus the argument that the cosmological time asymmetry may have origin in 
the geometric phase still survives. 
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